Conductance characteristics of current-carrying d-wave weak links 
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The local quasiparticle density of states in the current-carrying d-wave superconducting structures 
was studied theoretically. The density of states can be accessed through the conductance of the 
scanning tunnelling microscope. Two particular situations were considered: the current state of the 
homogeneous film and the weak link between two current-carrying d-wave superconductors. 
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INTRODUCTION 

Unconventional superconductors exhibit different fea- 
tures interesting both from the fundamental point of view 
and for possible applications In particular, double 
degenerated state can be realized in d-wave Josephson 
junctions If the misorientation angle between the 
banks of the junction x is taken tt/4, the energy minima 
of the system appear at the order parameter phase differ- 
ence </> = ±7r/2. These degenerate states correspond to 
the counter flowing currents along the junction bound- 
ary. Such characteristics make d-wave Josephson junc- 
tions interesting for applications, such as qubits ,3]. Our 
proposition was to make these qubits controllable with 
the externally injected along the boundary transport cur- 
rent [ij . It was shown that the transport current and the 
spontaneous one do not add up - more complicated in- 
terference of the condensate wave functions takes place. 
This is related to the phenomena, known as the param- 
agnetic Meissner effect [l|. 

It was demonstrated both experimentally [5( and the- 
oretically 0, Q that at the boundary of some high-T c 
superconductors placed in external magnetic field the 
current flows in the direction opposite to the diamag- 
netic Meissner supercurrent which screens the external 
magnetic field. This countercurrent is carried by the 
surface-induced quasiparticle states. These nonthermal 
quasiparticles appear because of the sign change of the 
order parameter along the reflected quasiparticle trajec- 
tory. Such a depairing mechanism is absent in the homo- 
geneous situation. Note that in a homogeneous conven- 
tional superconductor at zero temperature the quasipar- 
ticles appear only when the Landau criterion is violated, 
at v s > Aq/pp. Here v s is the superfluid velocity which 
parameterizes the current-carrying state, Ao stands for 
the bulk order parameter, and pf is the Fermi momen- 
tum. The appearance of the countercurrent can be un- 
derstood as the response of the weak link with negative 
self-inductance to the externally injected transport su- 
percurrent. The state of the junction in the absence of 
the transport supercurrent at zero temperature is unsta- 
ble at 4> = 7r from the point of view that small deviations 
5cf> = ±0 change the Josephson current from to its max- 



imal value The response of the Josephson junction 
to small transport supercurrent at <f> = it produces the 
countercurrent [9(. It is similar to the equilibrium state 
with the persistent current in ID normal metal ring with 
strong spin-orbit interaction: there is degeneracy at zero 
temperature and <f> = it, and the response of the ring is 
different at 6(f> ^ or B ^ 0, where B is the effective 
magnetic field which enters in the Hamiltonian through 
the Zeeman term (which breaks time-reversal symmetry) 
[ro| . The degeneracy is lifted by small effective magnetic 
field so that the persistent current rapidly changes from 
to its maximum value. In the case of the weak link 
between two superconductors in the absence of the trans- 
port supercurrent there is degeneracy between +p y and 
— p y zero-energy states; both the time-reversal symme- 
try breaking by the surface (interface) order-parameter 
and the Doppler shift (due to the transport supercurrent 
or magnetic field) lift the degeneracy and result in the 
surface (interface) current 

In recent years mesoscopic superconducting structures 
continue to attract attention because of the possible ap- 
plication as qubits, quantum detectors etc. (e.g. [3j, 
[111]). In particular, such structures can be controlled 
by the transport supercurrent and the magnetic flux 
(through the phase difference on Josephson contact). 
This was in the focus of many recent publications, e.g. 
H El, EI Q 0, (Hi- Here we continue to study the 
mesoscopic current-carrying cZ-wave structures. Particu- 
larly, we study the impact of the transport supercurrent 
on the density of states in both homogeneous film and in 
the film which contains a weak link. 



MODEL AND BASIC EQUATIONS 

We consider a perfect contact between two clean singlet 
superconductors. The external order parameter phase 
difference <f> is assumed to drop at the contact plane at 
x = 0. The homogeneous supercurrent flows in the banks 
of the contact along the y-axis, parallel to the boundary. 
The sample is assumed to be smaller than the London 
penetration depth so that the externally injected trans- 
port supercurrent can indeed be treated as homogeneous 



2 



far from the weak link. The size of the weak link is as- 
sumed to be smaller than the coherence length. Such 
a system can be quantitatively described by the Eilen- 
berger equation [8j. Taking transport supercurrent into 
account leads to the Doppler shift of the energy variable 
by p fV s . The standard procedure of matching the solu- 
tions of the bulk Eilenberger equations at the boundary 
gives the Matsubara Green's function G aJ (0) at the con- 
tact at x = (4J. Then for the component G* 1 = g(ui, r) 
of G w , which defines both the current density and the 
density of states (see below), we obtain in the left (L) 
and right (R) banks of the junction: 

9l,r (r) = g L ,R (oo)+[ff(0) - g L:R (oo)] exp f- -^j~r^ 

^ (1) 

gL,n(oo) = ——, (2) 



low temperature the dependence of the conductance on 
the bias voltage V is given by the following relation [13] : 

G(eV) = G N (D(p F )N{eV, p F )) , (8) 

where Gjy is the conductance in the normal state; D(pp) 
is the angle-dependent superconductor-insulator-normal 
metal barrier transmission probability. The barrier can 
be modelled e.g. as in Ref. Q with the uniform proba- 
bility within the acceptance cone < i? c , where $ is the 
polar angle and the small value of d c describes the thick 
tunnelling barrier: 

D(o)=i^-n (9) 

where #(...) is the theta function. 

CONDUCTANCE CHARACTERISTICS OF THE 
HOMOGENEOUS CURRENT-CARRYING FILM 



5(0) 



lo(£I l + fig) - isgn(^) A L A R sin < 
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(3) 



Here w = ttT(2ti + 1) are Matsubara frequencies, Al,a 
stands for the order parameter in the left (right) bank, 
and 
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(4) 



The direction-dependent Doppler shift p_pv s results in 
the modification of current-phase dependencies and in 
the appearance of the countercurrent along the boundary. 

The function g(u>,r) defines the current density, as fol- 
lowing: 



(5) 



w„>0 



Here iVo is the density of states at the Fermi level, (...)« 
denotes averaging over the directions of Fermi velocity 
vjr, v = vp/vp is the unit vector in the direction of vp. 
Analytic continuation of g(u>), i.e. 



ff(e) =ff(w -ie + 7), 



(6) 



gives the retarded Green's function, which defines the 
density of states: 



N(e,r) = Reg(e,r) 



(J) 



Here 7 is the relaxation rate in the excitation spectrum 
of the superconductor. 

The local density of states can be probed with the 
method of the tunnelling spectroscopy by measuring the 
tunnelling conductance G = dl/dV of the contact be- 
tween our superconducting structure and the normal 
metal scanning tunnelling microscope's (STM) tip. At 



Before studying the current-carrying weak link we con- 
sider the homogeneous situation. We will consider the 
d-wave film as shown in the left inset in Fig. [TJ The 
motivation behind this study is twofold: first, to demon- 
strate the application of the theory presented above, and 
second, to describe recent experimental results [14[. 
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FIG. 1: Normalized (divided by Gat) conductance dl/dV for 
the homogeneous current-carrying state in the <i-wave film 
for different values of the transport current. The curves are 
plotted with 7/A00 = 0.15 and § c = O.br [A o = A Q (v s = 0)]. 
Left and right insets show the schemes for probing the density 
of states in the current-carrying d-wave film and in the weak 
link (see text for details). 

The system considered consists of the <i-wave film, 
in which the current is injected along the y-axis, and 
the STM normal metal tip (another STM contact is not 
shown in the scheme for simplicity; for details see [3])- 
Following the experimental work [lj], we consider the 
c-axis along the x-axis and the misorientation angle be- 
tween a-axis and the direction of current (y-axis) to be 
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7r/4. Such problem can be described with the equations 
presented in the previous section as following [7l.ll5l]. 

Consider the specular reflection at the border, when 
the boundary between the current-carrying d-wave su- 
perconductor and the insulator can be modelled as the 
contact between two superconductors with the order pa- 
rameters given by = A($) = Aocos2(^ — \) an d 
A R = A(-tf) ee A and with = 0. Then from Eq. © 
we have the following: 



t3 



(n + o) 



AA' 



(10) 



where i7 = y/uj 2 + A 2 and £1 = \j uj 2 + A" . This expres- 
sion is valid for any relative angle \ between the a-axis 
and the normal to the boundary; in particular, 

g(u) = x = (A(0) = A cos2tf), (11) 

g(<j) = 2, X=7 (AO?)=A sin2tf). (12) 

The accurate dependence of the gap function A = 
Ao(w s , 7) can be obtained from Ref. [4| with introducing 
7 as following: p_fv s — ► p_fv s — 17 (which is analogous 
to Eq. ([6])). The energy values in this paper are made 
dimensionless with the zero-temperature gap at zero cur- 
rent: Aoo = A (w s = 0). 

And now with Eqs. $TZ§ and flUHJ) we plot the STM 
conductance for the current-carrying ci-wave film in Fig. 
(P) . We obtain the suppression of the zero-bias conduc- 
tance peak by the transport supercurrent, as was studied 
in much detail in Ref. [1J| . Our results are in agreement 



with their Fig. 1. Also the authors of Ref. [14| devel- 
oped the model based on phase fluctuations in the BTK 
formalism to explain the suppression of the zero-bias con- 
ductance peak. However, their theoretical result, Fig. 2, 
describes the experimental one only qualitatively, leaving 
several distinctions. They are the following: (i) position 
of the minima (eV/Aoo ~ 0.5 and 1 for the experiment 
and the theory respectively); (ii) height of the zero-bias 
peak at zero transport current (~ 2.5 and 4 respectively); 
(iii) height of the peak at maximal transport current 
(~ 1.3 and 2.5 respectively); (iv) presence/absence of 
the minima for all curves. Our calculations, Fig. (TTJ, 
demonstrate agreement with the experiment in all these 
features. The agreement we obtained with two fitting 
parameters, 7 and i? c . 

To further demonstrate the impact of the two fitting 
parameters of our model, 7 and i? C) in Fig. ([2|) we plot the 
normalized conductance fixing one of them and chang- 
ing another. The figure clearly demonstrates how they 
change the shape of the curves: the position of the min- 
ima, splitting of the zero-bias peak etc. Note that the 
splitting is suppressed at small $ c and high 7. This ab- 
sence of the splitting was observed in the experiment 14 1 
and studied in several articles, e.g. fljj ]. 
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FIG. 2: Normalized conductance dl /dV for the homogeneous 
current-carrying state in the d-wave film for different values 
of 7 and # e at pFV a / Aoo = 0.5. 



CONDUCTANCE CHARACTERISTICS OF THE 
CURRENT-CARRYING WEAK LINK 

Consider now the weak link between two d-wave 
current-carrying banks. For studying the effect of both 
the transport current and the phase difference on the 
density of states in the contact, we propose the scheme, 
presented in the right inset in Fig. [TJ The supercurrent 
is injected along y-axis in the superconducting film, as it 
was discussed in the previous section. Besides, the weak 
link is created by the impenetrable for electrons parti- 
tion at x — 0. The small break in this partition (a < £0) 
plays the role of the weak link in the form of the pinhole 
model @, [H The STM tip in the scheme is posi- 
tioned above the weak link to probe the density of states 
in it. Two more contacts along the x-axis provide the 
order parameter phase difference <f> along the weak link. 
This can be done, for example, by connecting the con- 
tacts with the inductance, as shown in the scheme, and 
applying magnetic flux $ e to this inductance. Then one 
obtains the phase control of the contact with the relation: 
<f> = $ e /$ - 

The two half-plains (for x < and x > 0) play the 
role of the two banks of the contact, which we also call 
left and right superconductors. In our scheme the banks 
carry the transport current along the boundary, and the 
Josephson current along the contact is created due to 
the phase difference. The banks we consider to be ci- 
wave superconductors with c-axis along the z-axis and 
with the misorientation angles \l = and \R = 7r /4- 
Now we can apply the equations presented in Sec. II to 
describe the conductance characteristics of the contact 
between current-carrying ci-wave superconductors. This 
is done in Fig. [3l where the normalized conductance is 
plotted for two values of the phase difference, for c/> = 
±7r/2 and with 7/A 00 = 0.1. The two values of the 
phase difference, <fi = are particularly interesting 

for the application since they correspond to the double- 
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degenerate states 0, 0). So, the density of states is the 
same in the absence of the transport supercurrent in both 
panels in Fig. [3] with mid-gap states (at eV < A 00 ) which 
create the spontaneous current along the boundary. The 
transport supercurrent (v a ^ 0) removes the degeneracy 
by significantly changing the mid-gap states (Fig. [3]), 
which explains different dependencies of the current in 
the contact on the applied transport current (i.e. on v s ), 
studied in 0,113]. 




FIG. 3: Normalized conductance dl/dV at the contact be- 
tween two current-carrying d-wave superconductors for dif- 
ferent values of the transport supercurrent (v 3 ) and for two 
values of the phase difference cj> = ±7i"/2. 



CONCLUSION 

We have studied the density of states in the current- 
carrying G?-wave structures. Namely, we have considered, 
first, the homogeneous situation and, second, the super- 
conducting film with the weak link. The former case was 
related to recent experimental work, while the latter is 
the proposition for the new one. The local density of 
states was assumed to be probed with the scanning tun- 
nelling microscope. The density of states at the weak link 
and the current (i.e. its components through the contact 
and along the contact plane) are controlled by the values 
of cf> and v s . The system is interesting because of possible 
applications: in the Josephson transistor with controlling 
parameters cf> and v s governed by external magnetic flux 
and the transport supercurrent [ll[, and in solid-state 
qubits, based on a contact of d-wave superconductors Q . 
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